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lO ■ Abstract 

The ACCMM model predicts the lepton spectrum from B meson 
(T) . decay by assuming the meson disintegrates into a spectator quark of 

definite mass and momentum distribution and an off shell b quark 

■ whose decay leptons (boosted into the rest frame of the meson) deter- 



mine the lepton spectrum. In this letter, we show that one can define a 



model dependent b quark mass so that the spectrum derived from the 
ACCMM model agrees very well with the free quark decay spectrum 
■ far from the endpoint. Near the endpoint, there is some disagreement, 

5_i ■ indicating the result is more model dependent. The integrated spec- 

tra are however very nearly identical. These results are in accordance 
with expectations based on the heavy quark effective theory. We con- 
clude that for LEP experiments, free quark decay might be as general 
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and more simple than the ACCMM ansatz for modeling the inclusive 
charged lepton spectrum from B meson decay. 
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1 Introduction 



Much attention has been devoted recently to the study of the lepton 
spectrum from inclusive semileptonic heavy hadron decay [|], ^, ||, [| . 
The key idea is that when a heavy hadron decays, the operator prod- 
uct expansion (OPE) and the heavy quark effective theory (HQET) 
may be employed to derive the decay spectrum. These methods were 
employed in the study of inclusive semileptonic decay Hi, — ► XIVl 
where H\, is a bottom hadron. The result (away from the endpoint of 
the spectrum) is that the inclusive differential decay rate ^ may be 
expanded in A/to, where A is a QCD related scale of order 500 MeV, 
and to is the mass of the heavy quark. The leading term (zeroth order 
in A/to) is the free quark decay spectrum, the subleading term van- 
ishes, and the subsubleading term involves parameters from the heavy 
quark theory, but should be rather small for b quark decay, as it is of 
order (A/ra) 2 . 

However, it is not always apparent how models reproduce this re- 
sult. For example, it is not manifest in the model of Altarelli, Cabbibo, 
Corbo, Maiani [||, and Martinelli, hereafter referred to as ACCMM, in 
which the differential decay spectrum of the charged lepton for inclu- 
sive semileptonic B meson decay is derived under some assumptions 
about the bound state. In this paper, we show that the linear terms 
in the differential distribution in the region far from the endpoint or 
of a suitably averaged spectrum do vanish (contrary to the claim in 
ref. Q). We present a detailed discussion of this model and show 
that by a suitable model dependent definition of the b quark mass, the 
linear terms in A/to may be eliminated so that the predictions of the 
ACCMM model agree well with those of the heavy quark effective the- 
ory (HQET). We first do the calculation directly from the ACCMM 
model, without invoking the formalism of HQET or the OPE. By ex- 
ploring this particular model, we hope to clarify the relation between 
models and the HQET result. Many of these ideas were discussed 
independently in ref. H . 
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Having shown that free quark decay reproduces the ACCMM spec- 
trum up to small corrections (of order (A/m) 2 which should be of the 
order of a percent), we discuss the question of whether simple free 
quark decay (with the b quark mass the free parameter) is as viable 
for modeling the lepton spectrum as the ACCMM model. We will 
see that the two models agree well far from the endpoint but deviate 
within a region of approximately 2pf from the endpoint. We consider 
the ACCMM spectrum in detail and also compare the results of an 
averaged ACCMM model to an averaged free quark model. We ask 
the question of whether smearing is required to get good agreement 
between the models, and if so, how much. It is useful to understand 
the endpoint region in a model in which there are no singularities 
introduced by using the OPE in the region where it does not con- 
verge. The examination of a particular model provides a complemen- 
tary approach to determining the necessary averaging which might be 
required in the endpoint region. Our conclusion is that free quark 
decay probably works as well for experiments at LEP, for example, as 
the ACCMM model and in fact has the same number of independent 
parameters (if the spectator mass of the ACCMM model is taken to be 
fixed). Whether the details of the endpoint region are important will 
however depend on the the application of the models and the accuracy 
which is required. 

Although we are considering free quark decay, there do exist higher 
order corrections. If accuracy better than a percent is required these 
higher order corrections must be included. These were considered 
in detail in refs. [^, EJ. The focus of this paper is not these small 
corrections, but to show that the linear correction terms can be made 
to vanish, and that free quark decay models fairly well the ACCMM 
spectrum (so long as the detailed structure of the endpoint is not 
important). However, we briefly comment in the conclusion on the 
generality of the higher order corrections of the ACCMM model. We 
also neglect perturbative QCD corrections in our discussion, but these 
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can be readily incorporated ||, [?[. Finally, we will consider a decay to 
the massive c quark, but our conclusions hold for a massless final state 
quark as well. The shape of the spectrum obtained in the massless 
case would however be significantly modified by perturbative QCD. 

We note that the ACCMM model was developed to consider in 
detail the endpoint of the lepton spectrum. As we will discuss, the 
ACCMM model is simply not presented in a form where its relation 
to the HQET result is obvious. In fact, it was explicitly constructed 
to avoid mention of a b quark mass. The LEP experiments appropri- 
ate this model, which was designed to study the endpoint, in order 
to estimate a systematic error in modeling the full spectrum. Our 
conclusion is that far from the endpoint, free quark decay is as good 
as the ACCMM model. The only model dependence comes in the 
treatment of the endpoint region. Any result depending on the details 
of this region is of course model dependent and unreliable. However, 
if it is only the distribution smeared over energy which is relevant, we 
will see that the ACCMM model and free quark decay give the same 
predictions at about the percent level. 

The outline of this paper is as follows. In the first section, we show 
how the ACCMM spectrum can be analytically expanded in b quark 
mass when far from the endpoint region. In the second section, we 
review the relevant results from the HQET and see how they apply to 
the ACCMM model. Section 4 considers in more detail the endpoint 
region of the spectrum, as well as the total integrated spectrum. We 
conclude in the final section. 

2 The ACCMM Model 

In order to incorporate the fact that when a b quark decays its spec- 
trum is not simply that from free quark decay, because the quark is in 
a bound state, the ACCMM model treats the decay of the B meson as 
a disintegration into a spectator quark of given mass and distribution 
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of momentum plus a b quark . The decay spectrum is determined by 
the kinematical constraints on the b quark. This should incorporate 
at least some of the corrections related to the fact that the b quark 
which decays is not free, but in a bound state. 

As we will discuss further in the following section, the heavy quark 
theory predicts that any model in accordance with its assumptions will 
give rise to a decay spectrum which up to small corrections (of order 
(A/m) 2 ) agrees with the free quark decay spectra, at least when suffi- 
ciently far from the endpoint (or when a suitable averaging procedure 
is applied). Since the ACCMM model does not obviously violate any 
of the underlying assumptions of a heavy quark model, we should ex- 
pect the resulting decay spectrum to very nearly agree with that of 
the decay of a free quark when far from the endpoint. Nonetheless, 
because the model does not incorporate the underlying QCD invari- 
ance, the predictions from this model are not compatible with the 
OPE/HQET result at the (A/m b ) 2 level. 

In ref. M, the results were obtained numerically, so it is difficult 
to identify the heavy quark mass dependence. Here, we do the heavy 
quark mass expansion explicitly. We first do the calculation in terms of 
the invariant mass of the b quark inside the B meson. We reexpress the 
answer in terms of the parameters used in the paper, namely ms, the 
B meson mass and pf, a parameter characterizing their assumed form 
for nonperturbative corrections. However, in this form, one cannot 
say anything about heavy quark mass corrections because we have not 
even defined the b quark mass. We show that if we define nib to be 
the average value of the energy of the b quark, the linear terms in the 
expansion of the differential distribution in powers of A/m do vanish. 
It is straightforward to generalize this result to arbitrary spectator 
and final state quark masses. 
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2.1 The Model 

The ACCMM model assumes a spectator quark of fixed mass m sp with 
a momentum distribution given by </>(|p|). For simplicity, this is taken 
to be Gaussian: 

which is normalized so that the integral over all momenta of 4>(,\p\)p 2 
(not only those which are allowed kinematically) is 1. 

Assuming conservation of energy and momentum tells us that when 
the momentum of the spectator quark has magnitude \p\ and mass 
m sp , the energy of the heavy quark will be 

E w = m B - \Jp 2 + m 2 sp (2) 

where tub is the B meson mass. 

The invariant mass of the b quark will then be 



W 2 = ml + m%-2m B ^ + ml (3) 

According to these assumptions, the distribution of lepton energy 
can be determined by boosting back the decay products of the b quarks 
of invariant mass W (in the b quark rest frame) to the rest frame of 
the meson and averaging over momenta. 

We define 

x = ^ (4) 



m 



(5) 



x m = 1 - e (6) 

Here, m/ is the final state quark mass. Then the lepton spectrum 
from the decaying b quark in its rest frame is 

dF(m b ,E) G 2 F mfx 2 (x m -x) 2 

~ltE— = 48vr3 (1-2)3 K 1 " ^ 3 " 2 *) + ( X " X ^ 3 ~ *)1 

(7) 
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Define 7 = Ew/W, (5 = p/Eyy- The spectrum of leptons in the 
ACCMM model is determined from 

dY B fPmax fld 2 T(W,E') . 

-§- = / dpp 2 4>{\p\) / ~ ,J , J dE' x 8 
dE Jo J 7 dE'd cos 

d COS On 



dcosO 



J ^pL 5{E - jE' - j(3E' cos 9 P ) (9) 



Here we have assumed the lepton is massless and we have defined the 
angle 6 p (associated with the distribution in momentum) with respect 
to the angle of the decaying lepton, for each of the orientations. It is 
straightforward to do the 9 P integral to derive 

d?B f r 2, /in 1 fd 2 T(W,E') 7 dE' . . 

We use the fact that (i^ 1 = pE\v/W 2 (p) and integrate over cos# to 
obtain 

dT B f, 2,n ix w2 [ Emax dE ' dT(W,E' 



dE J dpPm) 2p^L -W^dE^ (") 



where 



= ^ (12) 

£ mai = B — ( 1 - 7 yj I ( 13 ) 

2 \ {m B -m sp yj 

where the former are determined by the (5 function and the fact that 
the decay product comes from a quark of "mass" W and the latter 
comes from requiring that E max > E_ and is given for m sp = 0. 
This is the result of ACCMM (up to the 7 factor, which is a small 
correction). 
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2.2 Heavy Quark Expansion Applied to the 
Model 



We now consider the implications of the expression we have just de- 
rived. We first restrict our attention to massless final state and mass- 
less quarks. Furthermore, we consider the distribution far from the 
endpoint and assume small pf (so that E max = E + ) . We then have 

dT B rvm ax w 2 GlW 4 f E + dE' (2E\ 2 ( AE' 



dE 



= Jo m)p j E _ -W {wj I 3 - W) 



(15) 



We can explicitly evaluate the E' integral, to get 
dT B _ G\E 2 



(16) 

where it should be borne in mind that both W and Eyy are functions 
of p. 

We now consider this expression when pf is small. The value of 
Pmax can be taken to be approximately oo so far as the integral goes, 
making a negligible error. We also see that only p ^ mg will con- 
tribute significantly to the integral (in fact less than about 2pf), so 
that we can also expand the denominator in p/rriB- We expand in p 
and do the integrals, obtaining 

dr B _ dr% , G 2 F E 2 m 2 B (E \ Pf , „(p}\ dr. 



dE dE 127r 3 y / 7r \ m B J m B \m\) dE { ' 

where the first term represents the decay of a quark of mass m B . 

So if we expand in terms of the meson mass, it looks like the 
differential decay distribution is that for a free quark of mass m B 
which decays plus a linear correction term of order pf/m B . However, 
it is clear that in this form we learn nothing about the heavy quark 
expansion since the answer is expressed in terms of the meson rather 
than the quark mass. In fact, at this point, it is not even clear what 
the quark mass means in this model. Nonetheless, we can define the 
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mass of the b quark to be that mass for which this expression looks 
like free quark decay up to quadratic corrections. This definition is 
not as random as it sounds, since from equation ( |l6|) we see that the 
final spectrum is a function of W and Eyy. Therefore, the way to 
eliminate the linear terms in the decay distribution is to define the 
quark mass so that neither W nor E\y has linear corrections; that is 
define the quark mass m b by 

m b = (E w {p)) (18) 

Up to quadratic terms in p, we have 

(W) = (E w ) =m B - (p) = m b (19) 

Because 

(f(E w )) = f(m B ) - (p)f'(m B ) + 0(p 2 ) = f(m b ) + 0(p 2 /m 2 ) (20) 

we see that when we define the b quark mass this way that any func- 
tion f(Ew) which can be expanded in Pf/f(Eyy) will be equal to 
f(m b ) + 0(p 2 /m 2 ). It is easy to check that when we express the 
distribution as a function of m b there is no linear correction. From 
this viewpoint the vanishing of the linear term is not very deep. It 
is just the statement that all the p dependence is through Ew, so 
by defining m b to be its average, we eliminate the linear corrections 
to the free quark differential decay spectrum. Notice this is true even 
though the b quark is not on mass shell. In fact, we checked that if one 
constructed a nonrelativistic quark model with the b quark on shell 
and a Gaussian momentum distribution that the spectrum falls about 
halfway between the ACCMM prediction and a free quark model. 

This analysis makes clear why there is always a definition of b 
quark mass for which the linear correction terms vanish. Up to terms 
quadratic in p 2 /m 2 , the p dependence is all through E\y. Once we 
define the average of E\y to be the b quark mass, all linear corrections 
are eliminated (so long as we are in a regime where the expansion 
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in momentum is legitimate). Notice that in this model, the relation 
between the b quark mass and the quark mass defined at high energy 
can in principle differ by an amount of order Aqcd- However, this is 
an artifact of the model. 

It is easy to see that this analysis can readily be extended to the 
case when the spectator quark mass is nonzero. We then have 



(E w ) =m B - ^m 2 p + p 2 ) (21) 

Again, because everything can be expressed as a function of W (up 
to small corrections of order (A/m) 2 , there is only a single quantity 
we need to know. This is true despite the fact that naively, it appears 
that m S p and pf are independent parameters. We see that again, 
everything can be expressed in terms of {Ew)- Of course, at higher 
order in (A/m) 2 , this is no longer the case. This is briefly discussed 
in the conclusion. 

To extend the analysis to the case of nonzero final state quark 
mass, one can proceed as above and explicitly do the integral over E'. 
However, it is simpler to proceed as follows. We observe that 

E.-E. = ^ (22) 

which is small if p is small. Assuming that the integrand is a smooth 
function (ie we are far from the region near the endpoint where it 
varies rapidly) we have 

W 2 f E + dT q {W,E')dE' W 2 dT q (W,E) (E w +p 



r*+ dV q (W,E') dE' W dl q (W,E) 

Je^ dE' ~W ~ 2pE w dE ° 9 



2pE w J E _ dE' E' 2pE w dE a \E W - p 

l\ 2 dT q (W,p,E) 
1 ) dE {26) 

The p integral averages the differential distribution. Because it only 
depends on E\y (up to order p 2 /m 2 ), the result is the same as before; 
linear corrections vanish if we choose nib = {E\v)- Notice that in fact 
that any choice which differed from this one by terms of order p 2 /m 2 
would also suffice. (Here we chose the velocity of the b quark to be 
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v = (1,0,0,0).) The above argument is of course very general. As 
long as the function is varying smoothly and one is sufficiently far 
from the endpoint so that E max = E + , one expects corrections to the 
differential distribution to occur only at order (pf/nib) 2 - Again we see 
that it is model indpendent, that is independent of the detailed form 
of the momentum distribution, </>(|p|). 

Notice the same argument can be used to show that the ACCMM 
distribution agrees exactly with that for a free quark of mass mB—m sp 
for pf = 0. We therefore do not present results for pj = in what 
follows. 

We conclude that the differential distributions should agree well 
between the ACCMM model and free quark decay with the b quark 
mass determined as above. In Figs 1 and 2 we illustrate the agreement 
for pj = .15, .30, m sp = .15, mf ina i = 1.5. The discrepancy between 
the curves grows with pf. We choose these values for pj as they are 
the ones used in ||. 

Notice we could also have considered a massless final state quark. 
We have also checked our result for this case. Notice also that through- 
out this paper we are neglecting perturbative QCD corrections. We 
did however check that our integrated spectra agree with the numbers 
quoted in ACCMM when the QCD corrections are incorporated. For a 
massless quark, the detailed form of the QCD corrections needs to be 
incorporated. For a massive quark, it is essentially an overall factor. 

We will consider further the integrated spectra in Section 4. We 
first briefly review the HQET result. 

3 HQET Prediction for Spectrum 

The first paper to discuss the lepton spectrum in the context of the 
heavy quark theory is Chay et al. |[|]. By using the operator product 
expansion and then the heavy quark effective theory, they showed how 
to rigorously derive both the leading order QCD result (free quark 
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E/GeV 

Figure 1: The inclusive differential semileptonic decay rate (in units of 10 ^V^L) for the 
free quark model (solid curve) and the ACCMM model (dashed curve). The parameters are 
Pf = 0.15, m sp = 0.15, nif = 1.5. 
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Figure 2: The same as in fig. f for Pf = 0.3. 
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decay) and the correction terms. The correction terms are small (of 
order (A/m) 2 ) and depend on an unknown parameter, so our primary 
interest will be the leading and subleading order corrections. Most of 
this section is review from ref. [||, ||, Q . We include it to see when 
and why the ACCMM model should agree with free quark decay. 
Heavy quark decay proceeds via the Hamiltonian 

H = ^V a J?jl (24) 

where 

Jf = ^(l + 75 )6 (25) 

Jf = F,y(l +75 )r (26) 

where Vm is the KM matrix element . The dependence on the hadronic 
matrix element in the inclusive decay will be of the form 

W^ v = {2nfY J ^{PB-q-Px){B\j} 1 \X){B\J u \B) (27) 

i 

Now for fixed q 2 , one can study the following quantity 

T„M) = ~ij d'xe-^iBlTiJ^rmB) (28) 

The quantity has the property that the discontinuity across a 
cut gives W^y. That is, ImT^ = —ttW^ u . There is a cut in which 
extends between \[q 2 < q ■ v < (m 2 + q 2 — m 2 -,) / (2mb) where mc is the 
mass of the lightest hadron containing the c quark. There are other 
cuts corresponding to scattering processes which are irrelevant to this 
analysis. 

The amplitude T^ u can be perturbatively computed in a region 
far from singularities where the operator product expansion applies. 
The matrix element of the time ordered product of currents can be 
expanded in terms of the matrix elements of operators of the heavy 
quark theory. The coefficients of the operator are determined by eval- 
uating the matrix elements of the time ordered product between quark 
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and gluon states. The quantity itself is determined from the matrix 
elements of the heavy quark operators between the meson states. 

At any given order in AQCD/mb, T^ u is expanded in terms of a 
finite number of operators. The important conclusions which arise 
from this systematic expansion are 1) At order (A/m) , you repro- 
duce the free quark result, 2) At subleading order in (A/m), there 
are no corrections. This means the leading corrections to free quark 
decay arise only at order (A/m) 2 and are therefore small. In fact, in 
the heavy quark theory, there are only two possible operators for B 
mesons. However, one of the operators has an unknown matrix ele- 
ment, so the leading correction is not known, but has been evaluated 
in terms of the unknown parameter (see [|] ) . 

The key ingredient to the vanishing of the (A/m) correction was 
the use of the equations of motion. When these were used, there 
was an implicit assumption || that the heavy quark field, defined at 
leading order in (A/m) by 

b( x ) = e- im » v - x b v (x) (29) 

was defined in such a way that a potential mass counterterm van- 
ished. That is, the choice of m^ corresponds to the choice of coun- 
terterm (B\b v 5mb v \B) =0 0. There is a single physical parameter, 
namely the b quark mass. Notice that for spectator decays, every- 
thing depends on this quark mass (not the meson mass). This means 
for example that the rate is the same for the meson or the baryon 
containing a single b quark. 

However, there are limitations to the OPE approach. The OPE 
breaks down at the endpoint of the lepton spectrum. This is because 
the OPE does not converge at this point. This is apparent when 
higher order terms are included. Without smearing, the perturbative 
expansion in (A/m) does not work in the endpoint region. In terms of 
the dual picture, the endpoint region is the regime which one expects 
to be resonance dominated. That is, this is the regime where the 
spectrum depends on the details of the bound state. In order to obtain 
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predictions applicable to this region, a suitable smearing of energy 
must be applied. In either picture, one expects the range of energies 
which need to be averaged over to be of order Aqcd- In the next 
section, we explore these statements in the context of the ACCMM 
model. 

4 The ACCMM Spectrum 

In this section, we consider the spectrum of the ACCMM model in 
its entirety. We compare this model to the statements in the previous 
section. In particular, we will be interested to see when averaging is 
required, and if it is, how large a range of energy should be averaged to 
get good agreement? Of course the important question is really how 
well should we expect the predictions of different models to agree, 
given that the b quark is not really free. 

From the analysis of the previous section, we expect that for a 
sufficiently well averaged spectrum, we will get agreement with a free 
quark model (up to order (A/m) 2 ), but that the detailed forms of the 
spectrum near the endpoint will differ. Both these statements are true, 
as can be seen in Figures 1 and 2. In this section, we first consider the 
total integrated spectrum. We then discuss in detail how the ACCMM 
model differs from the free quark prediction. We averaged over energy 
to see how large a smearing is required to get good agreement (of 
course this depends on the accuracy desired). Finally, we discuss the 
question of how large an error in extracting physical parameters one 
is likely to make by using one model rather than the other. 

First let's consider the total integrated spectrum. From the anal- 
ysis of the previous section, we know that the integrated spectrum 
should agree well between the full result and our approximation of free 
b quark decay. This follows because the full integration is certainly 
averaged over energy over a sufficiently large interval. Therefore, the 
two should agree to within order (pf /tob) 2 . For example, for pj = .3, 
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m b 


X 


At 


A 


Pf = 0.1 


5.02 


1.733 


0.1 


0.5 


Pf = 0.2 


4.94 


1.929 


0.4 


2.1 


p f = 0.3 


4.84 


2.193 


0.8 


4.8 


Pf = 0.4 


4.73 


2.522 


1.2 


7.4 


p f = 0.5 


4.62 


2.926 


1.6 


10.86 



Table 1: For the values of pf in column 1, we display the corresponding value of J7l(, in the 
second column and the parameter X defined by T = X 10 14 V^BgL (f° r a quark of 
this mass, where BgL is the semileptonic branching fraction) in the third column. In the 
fourth column is given the fractional discrepancy in percent between the ACCMM determined 
lifetime and that from free quark decay, At = X 10 ^ V^^Bgj^. In the last column 
is given the deviation A = 100 fo' max I ^JJf ~ \dE/Y\. In all cases we have taken 

m sp = 0.15, mfinai = 1.5. 

the difference would be expected to be of order .4%. In Table 1, for 
various values of pf (with m sp fixed at .15), we give the corresponding 
values of m& and the lifetime according to the semileptonic decay of 
a quark of this mass. In the fourth column, we give the fractional 
deviation in percent , At between the ACCMM model determined 
lifetime and that of the free quark model. We see that there is very 
good agreement. 

It is easy to see that this had to be the case. Since the total 
spectrum for the ACCMM model is given by 

f dT B f , 9 , „ „ f dcos6L 1 

J dE' dF( ^f ) 6{E - -)'E' - 7 ftE' cos <)„) (30) 

which it is easy to see is the average of dT(W, E)/dE (up to quadratic 
corrections) . Again the average of this quantity agrees with the quan- 
tity evaluated at the average value up to higher order corrections. 
However, the entire spectrum does not agree so well. In particular, 
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within about 2pj of the endpoint, the spectra look different. In the 
fifth column is shown the ratio of the integrated absolute value of the 
difference of the differential spectra divided by the total integrated 
spectrum. We see that this discrepancy is relatively large. That is 
because the deviation can be of order Pf/m (note that it vanishes when 
Pf goes to zero where the distribution reduces to a delta function). 
It is only because the integrated spectrum averages the regions of 
positive and negative difference between the two models that the total 
integrated rates are in better agreement. 

We gain some insight into this cancellation by considering the exact 
expression for the differential decay rate, which is 

dT B f Pl , 2 ,n ^ W [ E+ dF(W,E') 



dE 



f Pl a 2, n ^ W f h+ dT(W,E') 

I d ^^Y P j E _ -V^ (31) 

rp max w fEi dT(W E') 

+ J iM -J E _ -U-L (32) 



9 

m ,. 



where E x = W(l - e)/2 and p x = m B /2(l f ) - E for m sp = 0. 

m 2 

This spectrum extends up to E max = mB 2 msp (1 - j^z^yi)- 

Notice that near the endpoint, (this is actually quite a large re- 
gion, of order at least 2pf), the result differs from the free quark decay 
spectrum due to two effects. First is that the term which would have 
reproduced the free quark decay spectrum (the first term) is no longer 
integrated over p up to p max ~ °o but is only integrated to p ma x — E. 
For E approaching the endpoint, the range of p integration is smaller 
than is necessary to get the full free quark decay spectrum. This re- 
duces the spectrum of the exact result relative to the free quark decay 
result as we approach the endpoint. However, as E approaches the 
endpoint, the second term gains significance. Because the spectrum 
is rapidly falling, the expansion in p followed by an average is not 
legitimate. Small values of p are favored, as these are equivalent to 
effectively larger quark mass. The integration over effective quark 
masses, some of which are larger than m^, permit the decay spectrum 
to continue beyond the endpoint of b quark decay. Near the endpoint, 
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the second effect increases the spectrum over the free quark decay 
result. 

Of course the spectrum near the endpoint is not well predicted, 
either by QCD or by models. It is for this reason that it is very 
unreliable. To extract a reliable value of m^, one could avoid this 
region and model the remaining measured parts of the spectrum with 
free b quark decay. Alternatively, one can average over energy. To see 
how large an energy interval is required (in this model), we averaged 
both the model and the free quark prediction using the averaging 
function 



Notice that an argument similar to that given above would show that 
so long as the averaged spectrum varies sufficiently slowly, one would 
expect the average (in momentum) of the smeared spectrum to be the 
smeared spectrum evaluted at the average W, so we should expect 
good agreement for sufficient smearing. 

The averaged spectra for pf = 0.3 and AE = 0.6 are displayed 
in Figure 3. We see that they agree rather well. This model yields 
some insight into how large a smearing region is necessary. In order 
to the endpoint to get results accurate at the level of about (pf/rrib) 2 , 
one must smear the result over a range of energies at least about 2pf. 
This can be seen in Table 2 where we give A for various values of 
the smearing parameter AE, for pf = 0.3. Because pf is a priori an 
unknown parameter, the required region is at least of order 500-1000 



However such precise agreement is perhaps not necessary. To 
understand how large an error is made by assuming a free quark 
model and neglecting higher order effects, we studied the question 
of how badly one would do in extracting a b quark mass assuming 
the ACCMM model were the correct description of a B meson. This 
gives some idea about how inaccurate the extraction of parameters 




(33) 



MeV. 
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Figure 3: The unaveraged inclusive differential semileptonic decay rates (in units of 10 ^fr) 
for the free quark and ACCMM model for Pf = 0.3, m sp = 0.15, TXlj = 1.5 and the 
smeared differential decay spectra for A.E = 0.6GeV. The upper two curves are the un- 
averaged ones;the solid curves correspond to the free quark model and the dashed ones to the 
ACCMM model. 





AE = 


AE = 0.3 


AE = 0.6 


AE = 0.9 


AE = 1.2 


A 


4.8 


2.7 


1.8 


1.4 


1.1 



Table 2: The deviation of the averaged spectra A (defined in Table 1) for different averaging 
values AE, in GeV. The fixed parameters are Pf = 0.3, m sp = 0.15, Ulf = 1.5. 
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-0.2 


-0.4 


-0.6 


-0.8 


AE = 


26 


20 


19 


22 


29 


AE = 0.3 


7.0 


3.9 


4.0 


7.3 


13 


AE = 0.6 


2.0 


0.85 


1.6 


4.4 


9.1 


AE = 0.9 


0.80 


0.27 


1.0 


3.0 


6.1 


AE = 1.2 


0.38 


0.12 


0.69 


2.0 


4.2 



Table 3: f° r the averaged and unaveraged spectra for Pj = 0.3 where \^ OC Jq max (/x — 
$2)^ f^dE, where j\ is a free quark spectrum and fi is the ACCMM spectrum. The different 
rows correspond to different averaging, and the different columns to the variation of the b 
quark mass (in percent). The absolute scale is arbitrary. 

would be by modeling the meson as free quark decay. We denned 
X 2 = $o max {fi ~ f2) 2 f2dEi, where f\ is a free quark spectrum and 
ji is the spectrum from the ACCMM model. We assumed the value 
of the quark mass one extracts would correspond to the minimum of 
this function, and varied the quark mass by fractional amounts given 
in percent in Table 3. The extracted b quark mass would always be a 
little smaller than the "correct" one, because the statistics are dom- 
inated by the region away from the endpoint where the quark model 
prediction is slightly less than that of the ACCMM model. Nonethe- 
less, the error in mass is quite small. In Table 3, we have taken 
Pf = 0.3 and smeared over different energies, AE. We see that even 
with no smearing, we would obtain a value of b quark mass which is 
accurate at about the 0.5% level. With smearing, one can do even 
better. With a smearing of AE = .6, one might obtain a b quark mass 
which is correct at the 0.2% level. Of course, there are no precise 
statistics behind this estimate. Nonetheless, it seems that free quark 
decay, even unaveraged, without higher order corrections, should de- 
scribe the b quark decay spectrum very well. Whether this is the case 
will of course depend on details and the required accuracy. 



21 



5 Conclusion 



A failing of the ACCMM model is that it does not incorporate prop- 
erly the parameters of QCD. If one fixes m sp (as was done in ref. ||), 
then the only real parameter of the model is pt. As we have empha- 
sized, this parameter is exchangeable for the b quark mass. But then 
there are no free parameters left to characterize Fermi motion and 
the spin dependent operator which occurs at second order. The sec- 
ond operator is not included at all; the first operator has a coefficient 
which is determined from pf, or equivalently m&. 

The relation between pj and the parameter K of ref [0, [| is 

w " <^f> " \% (34> 

From this point of view the values of K used in ref. [||, || correspond 
to fairly large values of pt . Of course the value of pj is unknown so 
these large values are potentially valid. The values of pf extracted in 
ref. H were determined not by the value of Fermi motion, but from 
the best quark mass to fit the data. In reality, Fermi motion requires 
an independent parameter. 

However, the ACCMM model does include one additional param- 
eter, namely m sp . Varying this parameter can change the relation 
between the b quark mass and pf. However, that the model does not 
incorporate the physics correctly is manifest at the level of (A/m) 2 
corrections. In the ACCMM model, the parameters (po/m 2 ) and 
((p) 2 /m 2 ) occur as independent parameters, both giving rise to or- 
der (A/m) 2 corrections. However, in a gauge invariant formulation, 
the first matrix element would correspond to the expectation value of 
Dq which is mass suppressed by the equations of motion. Gauge in- 
variance imposes restrictions which are not necessarily true in a model. 
This will presumably be a problem of any model which does not in- 
corporate gauge invar iance and in which (j?q) does not vanish. 

Nonetheless, we see from Table 4, that independent of the value 
of m S p, the free quark decay model agrees well with ACCMM (in fact 
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m h 


X 


At 


A 


m sp = 0.15 


4.84 


2.193 


0.82 


4.84 


m sp = 0.3 


4.75 


2.462 


1.21 


4.34 


m sp = 0.5 


4.60 


3.040 


1.67 


4.36 



Table 4: The same as in table 1, for different values of Tfl S p. The fixed parameters are 
Pf = 0.3, rrif = 1.5. 

the deviation from ACCMM is not strongly dependent on m sp ). This 
is not surprising, since for small m sp , it is obviously not important, 
while for large m sp , it looks more like a nonrelativistic quark model. 

In summary, we have shown that the b quark mass can be defined 
so that the predictions of the ACCMM model agree with that of a free 
quark model at the level of (A/to) 2 , so long as one is sufficiently far 
from the endpoint or a suitable energy average is applied. Only for 
results which depend on the detailed form of the endpoint region would 
this approximation be inadequate. However such results would be 
unreliable in any case. For the purpose of LEP experiments, free quark 
modeling might be sufficient . The data is presumably effectively 
smeared by the fragmentation function of the B meson, even without 
explicit smearing. 

In our analysis, we have ignored perturbative QCD corrections 
and higher order effects. The first are straightforward to incorporate 
[||, ||. To incorporate the latter as generally as possible would require 
leaving the Fermi motion parameter as independent, and including 
the known spin dependent operator. This would allow one to include 
higher order effects if greater precision is necessary. The error in 
using a free quark model can presumably be estimated by varying the 
amount of smearing and the unknown higher order parameter. 
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